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ABSTRACT 

We present a case study for the relevance of the r-mode instability for accreting compact 
stars. Our estimates are based on approximations that facilitate back-of-the-envelope calcula- 
tions. We discuss two different cases: 1) For recycled miUisecond pulsars we argue that the 
r-mode instability may be active at rotation periods longer than the Kepler period (which pro- 
vides the dynamical limit on rotation) as long as the core temperature is larger than about 
2 X 10^ K. Our estimates suggest that the instability may have played a role in the evolution 
of the fastest spinning pulsars, and that it may be presently active in the recently discovered 
2.49 ms X-ray pulsar SAX J1808. 4-3658 as well as the rapidly spinning neutron stars observed 
in low mass X-ray binaries. This provides a new explanation for the remarkably similar rotation 
periods inferred from kHz quasi-periodic oscillations in the LMXBs. The possibility that the 
rotation of recycled pulsars may be gravitational-radiation limited is interesting because the 
gravitational waves from a neutron star rotating at the instability limit may well be detectable 
with the new generation of interferometric detectors. 2) We also consider white dwarfs, and 
find that the r-mode instability may possibly be active in short period white dwarfs. Our order 
of magnitude estimates (for a white dwarf of M = Mq and R = O.Oli?© composed of C^^) show 
that the instability could be operating for rotational periods shorter than P ^ 27 — 33 s. This 
number is in interesting agreement with the observed periods (> 28 s) of the rapidly spinning 
DQ Herculis stars. However, we find that the instability grows too slowly to affect the rotation 
of these stars significantly . 
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1. Introduction 



dwarfs. 



The recent discovery of a new instability in rotating rel- 
ativistic stars has revived interest in problems concerning 
mechanisms that may limit the rotation rate of neutron 
stars. As was first shown by one of us (Andersson 1998), 
the so-called r-modes (Papaloizou and Pringle 1978) are 
formally unstable at all rates of rotation in a perfect fluid 
neutron star (see also Friedman and Morsink 1998). The 
mechanism behind the instability is the familiar one that 
was discovered by Chandrasekhar (1970), Friedman and 
Schutz (1978): The modes are unstable due to the emis- 
sion of gravitational waves (see Stergioulas (1998) for a 
review on nonaxisymmetric instabilities in rotating rel- 
ativistic stars). Even more surprising than the actual 
existence of the instability are the results of recent es- 
timates of the time-scale at which the r-mode instability 
will grow (Lindblom, Owen and Morsink 1998; Andersson, 
Kokkotas and Schutz 1998). It seems that the r-mode in- 
stability has the potential to slow down newly born neu- 
tron stars dramatically, and that the r-mode instability is 
considerably stronger than the previously considered one 
for the f-mode of the star (cf. the results of Lindblom 
1995; Stergioulas and Friedman 1998). The new results 
suggest that, assuming that a neutron star is born spin- 
ning as fast as it possibly can (at the Kepler frequency 
~ O.Q7\/ttGp, i.e. with a period Pk ~ 0.5 — 2 ms), 
the r-mode instability will force the star to spin down to 
a period of roughly 15 ms in the first year or so. This pre- 
diction is in good agreement with the inferred initial pe- 
riod for the Crab pulsar (19 ms). The uncertainties in the 
present models may also accomodate the newly discovered 
16 ms pulsar in N157B (Marshall et al 1998). The birth 
spin period for this rapidly spinning pulsar is less certain, 
but one can estimate that it ought to have been shorter 
than 9 ms (this follows by assuming the braking index to 
be in the range observed for young pulsars together with 
the youngest likely age of 5000 years). Anyway, if the 
r-mode instability is active in a newly born neutron star 
it will force the star to spin down. The large fraction of 
the initial rotational energy that is then radiated away 
as gravitational waves may well be detectable with, for 
example, LIGO (Owen et al 1998). 

These are undoubtedly exciting suggestions, that indi- 
cate that the r-mode instability plays an important role in 
astrophysics. Given the promise of the suggested scenario 
for newly born neutron stars, we are inspired to investi- 
gate whether the instability can be relevant also in other 
situations. In this paper we assess the potential relevance 
of the instability for accreting stars, both neutron stars 
and white dwarfs. Specifically, we provide rough estimates 
of the timescales involved and ask whether the r-mode in- 
stability provides an upper limit on the rotation of stars 
spun up by accretion. We compare our estimates to ob- 
served data for the recycled millisecond pulsars (MSPs) 
and the rapidly spinning neutron stars in low-mass X- 
ray binaries (LMXBs), as well as the DQ Herculis white 



2. Back of the envelope timescales 

At the present time there are no fully relativistic stud- 
ies of the dynamics of perturbed rotating stars (although 
the neutral modes that signal the onset of instability 
of the f-mode have been calculated, see Stergioulas and 
Friedman (1998) for a fully relativistic calculation and 
Yoshida and Eriguchi (1997) for a calculation in the rel- 
ativistic Cowling approximation). Until such calculations 
become available we will not know the detailed effects of 
a gravitational-wave instability. However, we can obtain 
useful estimates that should be qualitatively correct, and 
may not differ too much from the quantitative truth, in a 
straightforward way (Lindblom et al 1998; Andersson et 
al 1998). In this paper we will use the estimates obtained 
by Kokkotas and Stergioulas (1998) in the context of uni- 
form density stars to assess the relevance of the r-mode 
in various situations. We prefer to work within this ap- 
proximation since the way that the various results scale 
with the parameters of the star (mass and radius) is then 
clear. 

The uniform density approximation leads to the grav- 
itational radiation growth time: 

_ /1.4AfQ\ /10km\^ f P V 

for the I — m = 2 r-mode (that leads to the strongest 
instability). M, R and P are the mass, the radius and 
the period of the star, respectively. This estimate of tg^ 
is roughly a factor of two smaller than the results for the 
iV = 1.0 polytropic equation of state and the specific stel- 
lar parameters {M = 1.5Mq and i? w 12.5 km) considered 
by both Lindblom et al (1998) and Andersson et al (1998). 
In a way, this discrepancy illustrates the uncertainties as- 
sociated with the present studies of the r-mode instability. 
Furthermore, since models of realistic neutron stars typi- 
cally have a mean polytropic index in the range 0.5 — 0.8, 
the difference of a factor of two or so (between polytropes 
of index and 1) in the gravitational radiation timescale 
can be seen to represent the uncertainty in the real equa- 
tion of state of neutron star matter. Also, as we will argue 
below, the difference of a factor of two has little effect on 
the conclusions of the present discussion. Consequently, 
we will work with the uniform density approximation (|l|) 
in the remainder of this paper. 

The damping timescale due to shear viscosity in a nor- 
mal fluid star is 

^--^■2^10 (^j (lOl^j [Wk) ^' 

(2) 

where T is the temperature of the star, and we have used 
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the shear viscosity coefficient 

?7„/ = 2x10^2 f (t^^ gcni^^s^^ 

^ VlQiSg/cm^ ViO'K; 

(3) 

This assumes that the main contribution to the shear 
viscosity arises from neutron-neutron scattering, which 
should be a vahd approximation at temperatures above 
that which the star becomes superfluid (Flowers and 
Itoh 1979). In a similar way to the gravitational- wave 
timescale, the normal fluid shear viscosity timescale (H) 
has a factor of two uncertainty. 

When the outer core of the star has become superfluid 
the main contribution to the viscous dissipation is due to 
electron-electron scattering. Then the relevant viscosity 
coefficient and timescale are: 

Vsf = 6 X 1022 ( P -\ ( ' g cm-1 s-i , 

(4) 

and 

For a superfluid star we should also be concerned with 
the so-called mutual friction. Specifically, it has been sug- 
gested that mutual friction will completely suppress the 
instability of the f-mode (Lindblom and Mendell 1995). 
At present there is no similar calculation for the r-mode 
instability, but it would not be surprising if mutual friction 
also has a strong effect on the r-modes. However, Mendell 
(1991) has suggested that mutual friction is dominated by 
the shear viscosity [as given by (||)] at temperatures below 
about 10^ K. This would indicate that the effect may not 
be dominant in the oldest stars. Given that results for the 
effect of mutual friction on the unstable r-modes are out- 
standing we will not include this effect in the present dis- 
cussion. For superfluid stars we will only account for the 
shear viscosity due to electron-electron scattering. This 
means that our results may overestimate the strength of 
the instability in the superfluid case. 

On the other hand, one would not expect the entire 
star to become superfluid: A cold neutron star will have 
an extended solid crust. This may be highly relevant for 
the r-mode estimates: The mode used in the calculation 
of the estimated timescales tg^ and tgv is in many ways 
similar to the f-mode of the star (albeit with a toroidal 
rather than spheroidal angular dependence). Specifically, 
both modes are well described by eigenfunctions that grow 
as r' towards the surface (in this sense they are "global" 
modes). For this reason, we expect them to be affected in 
a similar way by the presence of a crust, and as was shown 
by McDermott et al (1988) the crust has hardly any effect 
on the frequency of the f-mode (basically since the mode 
propagates at the speed of sound which is generally much 
higher than the speed of a shear wave in the crust). We 



expect that this is likely to hold also for the r-mode under 
consideration, and therefore the above timecales may, in 
fact, be reasonable also for a star with a superfluid core 
and a sohd crust. 

As long as we are mainly interested in relatively low 
temperatures, and other, more exotic, dissipation mecha- 
nisms do not play a role, the effect of the r-mode instabil- 
ity can be inferred from tgw and tsv ■ At high temperatures 
(above 10^ K) one expects the bulk viscosity of the fluid 
to play a significant role, cf. the estimates of Lindblom et 
al (1998) and Andersson et al (1998). However, in accret- 
ing systems the temperature should always be low enough 
that it is sufficient to include only tg^ and . When this 
is the case one can easily deduce the critical rotation pe- 
riod above which the r-mode instability affects the rota- 
tion of the star (when the gravitational- wave growth time 
is shorter than the viscosity damping time) . We thus find 
that the mode grows if the period is shorter than 

(6) 

for a normal fluid star, and 

/ D \ 3/2 / rp \ 1/3 

when we use the viscosity due to electron-electron scatter- 
ing in a superfluid. Interestingly, these critical periods are 
not strongly dependent on the mass of the star. Further- 
more, the uncertainties in (|l]) , (^) and (||) have little effect 
on the critical period. For example, the uncertain factors 
of two in tgw and tsv individually lead to an uncertainty 
of 12% in Pc- When combined, the uncertainties suggest 
that we may be (over)estimating the critical period at the 
25 % level. Considering uncertainties associated with the 
various realistic equations of state for supranuclear mat- 
ter, and the many approximations upon which our present 
understanding of the r-mode instability is based, we feel 
that it is acceptable to work at this level of accuracy. 

3. Implications for millisecond pulsars 

We will now discuss the possibility that the r-mode in- 
stability may be relevant for the period evolution of the 
fastest observed pulsars. All observed millisecond pulsars 
have periods larger than the 1.56 ms of PSR1937-I-21, and 
it is relevant to ask whether there is a mechanism that pre- 
vents a neutron star from being spun up further (eg to the 
Kepler limit) by accretion. Specifically, we are interested 
in the possibility that the r-mode instability plays such a 
role. Before proceeding with our discussion, we recall that 
Andersson et al (1998) have already pointed out that the 
instability has implications for the formation of millisec- 
ond pulsars (albeit in an indirect way). Specifically, the 
strength of the r-mode instability seems to rule out the 
scenario in which millisecond pulsars (with P < 5 — 10 ms) 
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are formed as an immediate result of accretion induced 
collapse of white dwarfs. Continued accretion would be 
needed to reach the shortest observed periods. In other 
words, all millisecond pulsars with periods shorter than 
(say) 5-10 ms should be recycled. 

Our main question here is whether it is realistic to ex- 
pect the instability to be relevant also for older (and in 
consequence much colder) neutron stars. Even though the 
critical period is much shorter for a cold star our estimates 
(H) and (^ are still above the Kepler period (« 0.8 ms 
for our canonical star) which suggests that the instability 
could be relevant. As an attempt to answer the question 
we will confront our rough approximations with observed 
data for MSPs, and the neutron stars in LMXBs. 

3.1. The millisecond pulsars 

In this section we discuss the r-mode instability in the 
context of the recycled MSPs. These stars are no longer 
accreting, and supposing that they have been cooling for 
some time they should not be affected by the instability at 
present. Our main question is whether the observed data 
is in conflict with a picture in which the r-mode instability 
halted accretion driven spin-up at some point in the past. 

First of all, our estimates show that the rotation will be 
limited by the Kepler frequency (using Pk ~ 0.8 ms for 
a canonical star) if the interior of the star is colder than 
T « 2 X 10^ K. Also, it is straightforward to show that, 
in order to "rule out" the instability (to lead to a critical 
period equal to the Kepler period at e.g. temperature 
4 X 10^ K) the dissipation coefficient of the shear viscosity 
(or any other dissipation mechanism) must be almost six 
orders of magnitude stronger than (|^) . 

Our inferred critical periods, (||) and (^ (for a canon- 
ical neutron star) are illustrated, and compared to ob- 
served periods and upper limits on the surface temper- 
atures (from ROSAT observations; see data given by 
Reisenegger (1997)) for the fastest millisecond pulsars in 
Figure In the figure we also indicate the associated 
upper limits on the core temperatures, as estimated using 
eqn (8) of Gudmundsson, Pethick and Epstein (1982). 

The illustrated r-mode instability estimates would be 
in conflict with the MSP observations if the interior tem- 
perature of a certain star were such that it was placed 
considerably below the critical period for the relevant tem- 
perature. Basically, an accreting star whose spin is limited 
by the r-mode instability would not be able to spin up far 
beyond the critical period, since the instability would ra- 
diate away any excess accreted angular momentum. As 
the accretion phase ends, the star will both cool down and 
spin down (the timescales for these two processes, photon 
cooling and magnetic dipole braking, are such that a MSP 
would evolve almost horisontally towards the left in Fig- 
ure |l|). 

Given the uncertainties in the available data we do not 
think the possibility that the r-mode instability may have 



played a role in the period evolution of the fastest MSPs 
can be ruled out. First of all, it must be remembered that 
the ROSAT data only provide upper limits on the surface 
temperature, and the true temperature may well be con- 
siderably lower than this. If the true core temperatures 
of the fastest spinning pulsars were roughly one order of 
magnitude lower than the ones indicated in Figure |l| our 
inferred lower limit on the spin period would, in fact, be in 
quite good agreement with observations. One must also 
recall the many simplifying assumptions that went into 
our estimates of the critical period. (Recall that we have 
already suggested that the uncertainties in the critical pe- 
riod may be at the level of 25%.) The most crucial of these 
assumptions regards the detailed role of superfluidity. It 
is not at all unlikely that superfluid mutual friction, or 
some other exotic dissipation mechanism, serves to sup- 
press the instability and decrease the critical period. In 
conclusion. Figure |l| is suggestive. It seems plausible that 
the r-mode instability is relevant for neutron stars spun 
up to millisecond periods. Furthermore, the figure indi- 
cates how improved temperature observations should be 
able to test present and future (more detailed) models for 
the r-mode instability. 

3.2. Neutron stars in LMXBs 

In the standard scenario it is assumed that MSPs have 
evolved from LMXBs, in which they spin up to the present 
rotation rate by accreting matter from their low-mass 
(< 0.4Mq) companions. If the accretion rate is close to 
the Eddington limit (lO^^A/Q/yr) it would typically take 
10^ years to spin a neutron star up to a period of a few 
milliseconds. Recent observations have provided interest- 
ing support for this model. First of all, the observed kHz 
quasiperiodic oscillations (QPOs) detected by the Rossi 
X-ray timing explorer (see van der Klis (1997) for a re- 
view) show that the neutron stars in LMXBs are spinning 
rapidly. Secondly, the recent discovery of the 2.49 ms X- 
ray pulsar SAX J1808. 4-3658 (Wijnands and van der Klis 
1998) provides the first evidence of an evolutionary link 
between the LMXBs and the MSPs. 

The observations of kHz QPOs from systems contain- 
ing rotating neutron stars provide a wealth of interesting 
information. Of specific relevance for our present discus- 
sion is the suggestion that many neutron stars in LMXBs 
are spinning with almost identical periods. This sugges- 
tion follows from the fact that the separation between the 
two KHz QPOs that are commonly observed stays con- 
stant even though the individual peak frequencies vary. 
The proposed explanation — the beat-frequency model 
(see van der Klis (1997) for references) — is that this ro- 
bust frequency is the underlying rotation frequency (or 
a multiple thereof). Assuming that this is the case one 
finds that the observed QPOs indicate a typical rota- 
tion frequency in the range 260-330 Hz {P in range 3- 
3.8 ms). (Interestingly, Backer (1998) recently suggested 
that MSPs are typically born with a period of 3 ms af- 
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ter spin-up.) There is, however, one major caveat to this 
interpretation of the QPOs: For Sco X-1 the peak separa- 
tion varies considerably with the luminosity of the system. 
It has, in fact, been argued that the available data is con- 
sistent with such a variation in all LMXBs (Psaltis et al 
1998). If this is the case the peak separation cannot corre- 
spond to the rotation frequency (since the rotation of Sco 
X-1 would then have to vary between 3.2 ms and 4.4 ms 
on short timescales) , but even so it seems reasonable that 
the observed beat-frequency is close (but not equal) to the 
rotation frequency. Hence, we believe it makes sense to 
assume that the neutron stars in the LMXBs are actually 
spinning at almost identical periods and try to understand 
why this is the case. 

Two explanations for this phenomenon have been pro- 
posed. In the first model (White and Zhang 1997) the star 
reaches an equilibrium period when the spin-up torque 
due to accretion is balanced by the spin-down torque due 
to the magnetic field. For a magnetic field dominated by 
the dipole the standard equilibrium period follows from 

Within this model one can use the observed periods and 
luminosities (which imply the accretion rates) and deduce 
a maximum magnetic field for each star. This results (cf. 
Table 1 of White and Zhang (1997)) in magnetic fields 
very similar to those found for MSPs. This is suggestive, 
but the magnetic fields deduced for the LMXBs will only 
compare favourably to the observed values for the MSPs if 
the field does not decay as the neutron star evolves. How- 
ever, this is in accord with the favoured model at present: 
The magnetic field is only expected to decay considerably 
during the accretion phase (Romani 1990; Urpin, Geppert 
and Konenkov 1998). Hence, the values that White and 
Zhang (1997) infer for B are reasonable as long as the low- 
field systems are close to the end of the accretion phase 
(which could well be the case since they are the ones with 
the lowest X-ray luminosity). 

However, the model of White and Zhang only leads to 
similar equilibrium periods if there is some intrinsic asso- 
ciation between M and B"^ [see (^)] and, as was pointed 
out by Bildsten (1998), no such relation is seen in the 
data for the (slower rotating) highly magnetized X-ray 
pulsars. In view of this, Bildsten proposed an alterna- 
tive explanation for the similar rotation periods in the 
LMXBs. He argued that the rotation may be limited by 
gravitational radiation in these systems. Specifically, he 
suggested that temperature gradients in the star lead to 
a time-dependent quadrupole moment, i.e. emission of 
gravitational waves that balance the spin-up torque due 
to continued accretion. This is an interesting suggestion, 
especially since the resultant gravitational waves may well 
be detectable. 



It is interesting to compare the data for the LMXBs to 
our inferred critical period for the r-mode instability. In 
order to do this, we need to estimate the interior tempera- 
ture of the neutron stars in LMXBs. For a star that is ac- 
creting close to the Eddington limit hydrogen (and possi- 
bly helium) burning of the accreted material can heat the 
core of the star up to a temperature considerably higher 
than the typical one for non-accreting neutron stars (Fu- 
jimoto et al 1984; Brown and Bildsten 1998). The sug- 
gested core temperature for a rapidly accreting neutron 
star is in the range 1 — 3 x 10* K. We use this estimate 
as an upper limit and compare the resultant data to our 
critical periods in Figure |l|. The Fi gure suggests that the 
r-mode instability provides a possible explanation for the 
inferred periods in the LMXBs. 

Regarding the comparison between our critical periods 
and the LMXB data in Figure |^ it is worthwhile making a 
few remarks. First of all, it is not impossible that the in- 
terior temperature of the LMXB neutron stars is actually 
lower than the 1 — 3 x 10* K we used as our upper limit. 
As was shown by, for example, Fujimoto et al (1984) a 
pion condensate in the core of the star can act as a heat 
sink and keep the bulk of the star at 10^ K or so, even 
though the outer layers of the star are heated up further 
by hydrogen burning. It is relevant to point out that our 
critical period would be in perfect agreement with such a 
model. Secondly, it is worth pointing out that the shear 
viscosity must be 200 times stronger than our superfluid 
value (^ in order for the r-mode instability not to be ac- 
tive in these systems. If this were the case we would have 
Pc ~ 3 ms at r « 3 X 10* and the r-modes would be stable 
in all indicated LMXBs. Finally, it should be noted that 
the r-mode explanation is in qualitative agreement with 
a possible trend seen in the LMXB: The neutron star can 
only spin up further as the accretion slows down. In our 
picture (and also in Bildsten's explanation) the star must 
accrete slower, and cool down, in order to be able to spin 
up. This means that the evolution towards periods shorter 
than 3 ms should take place on a cooling time-scale. In 
contrast, in the model of White and Zhang (1997) the star 
would spin up on the timescale that the magnetic field de- 
cays. In reality, these two timescales are, of course, likely 
to be linked in some way. 

Having suggested an alternative explanation for the 
clustering of the observed rotation periods in LMXBs 
(that the r-mode instability prevents the neutron stars 
from spinning up once they reach the critical period) one 
would perhaps expect us to attempt to rule out the other 
two models. However, given the uncertainties in all the 
suggested scenarios, we do not believe it is meaningful to 
rule out either possibility at the present time. To be able 
to do so we need more detailed modelling. In fact, we can 
see no reason why several of these mechanisms could not 
be operating at the same time, or be relevant in differ- 
ent individual systems. As one can immediately deduce 
from (0) and (||) the final outcome depends on a delicate 
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balance. For example, in order for the r-mode instability 
to dominate over the magnetic dipole we need Peq < Pc- 
Using (0) and (^) we can translate this into a relation be- 
tween the temperature, the accretion rate and the mag- 
netic field strength. From this relation we find (for our 
canonical star): i) For accretion at the Eddington rate and 
a temperature of 10^ K, the r-mode instability dominates 
for B < 3 X 10^ G, i.e. for most of the observed systems, 
ii) For a lower accretion rate, 1% of the Eddington limit, 
at 10^ K we find that the r- modes would be dominant for 
< 1.3 X 10* G, i.e. only a few of the observed systems. 
This simple example shows that the two mechanisms are 
active at roughly the same level and that neither of them 
should be ruled out at the present time. 

3.3. Gravitational-wave estimates 

As we have seen, the rotation of the fastest observed 
neutron stars in accreting systems may be limited by the 
r-mode instability. Hence, it is meaningful to estimate 
the amplitude of the gravitational waves that could be 
emitted from these systems. As was first pointed out by 
Wagoner (1984), a gravitational-wave instability in an ac- 
creting star will lead to periodic gravitational waves (as- 
suming that the system reaches a steady state where the 
excess angular momentum of the accreted matter is radi- 
ated as gravitational waves). Following Wagoner (1984) 
we can deduce the gravitational-wave amplitude by using 
the angular momentum gained by accretion: 



j ^ VGMRM = E , 



OJ 



(9) 



(where uj is the angular frequency of the mode) in the 
standard flux formula 



4G / 1 



\E\ 



(10) 



where r is the distance to the source, and h is the (di- 
mensionless) gravitational-wave amplitude. Using these 
relations, together with our r-mode estimates, we find 



2.3 X 10 




To illustrate what this means, let us provide a spe- 
cific example. Since it would provide one of the strongest 
sources we consider the particular case of Sco X-1 (cf. 
Bildsten 1998). Assume that the true rotation period of 
Sco X-1 is 4 ms, and that the interior temperature is such 
that the r-mode instability is active and provides a limit 
on the spin of the star (this corresponds to assuming that 
T « 10* K cf. Figure]^). Given this, and the relevant data 
for Sco X-1 (r = 0.7 kpc and an average accretion rate of 
M « 3 X IQ-^Mq/yi (Lang 1992)) we find that (0) im- 
plies h « 3.5 X 10^^® for a star with canonical mass and 



radius. The effective amplitude achievable after matched 
filtering improves as the square-root of the number of de- 
tected cycles, so we get /ig// — ^/ujt^^b^j2nh ^ 10~^^ 
after about two weeks of integration. After one year of 
integration, a similar source at the centre of our galaxy 
(at 10 kpc) would reach the same effective amplitude. In 
other words, the gravitational waves from galactic neu- 
tron stars whose rotation is limited by the r-mode insta- 
bility should be observable by, for example, LIGO. The 
detection of this kind of essentially continuous signals is 
obviously more complicated than we have indicated here, 
but since the data analysis issues are identical to those 
for gravitational waves from a slightly deformed rotat- 
ing neutron star we refer the interested reader to detailed 
studies of that problem (Jaranowski, Krolak and Schutz 
1998, Brady et al 1998). 

Interestingly, given a detection of gravitational waves 
from the LMXBs one would easily be able to distinguish 
between the r-mode instability and radiation from the 
time-dependent quadrupole model suggested by Bildsten 
(1998). The r-modes radiate at 4/3 of the rotation fre- 
quency, while a rotating quadrupole should generate ra- 
diation of (mainly) twice the rotation frequency. Suppos- 
ing that the rotation frequency is known from a QPO one 
should have no trouble distinguishing between the two 
mechanisms. 

Finally, we should also mention that there may exist 
systems that are radiating gravitational waves according 
to the prescribed scenario, but are otherwise invisible. As 
was suggested by Schutz (1997) the accretion rate may 
be close to the Eddington limit if the neutron star is part 
of a so-called Thorne-Zytkow object, i.e. when it spirals 
inside a giant star at the late stages of binary evolution. 
Such systems should be similar to the example we dis- 
cussed above, and would likely be detectable through the 
radiated gravitational waves. 

4. Implications for white dwarfs 

Although it is clear that the r-modes would formally be 
unstable also in less compact perfect fluid stars, it is not 
at all obvious that the timescales will work out in such 
a way that the instability has any relevance. Hence, it 
is interesting to attempt to extend our estimates to white 
dwarfs. To achieve a rough estimate for the critical period 
is rather straightforward. We first rescale our expression 
for the growth time-scale due to gravitational waves (|l|) 
to canonical values for a white dwarf. This leads to 



3.1 X 10-' 



M 



OMR 



■0 



R 



— 1 yi-s . (12) 

30s' ^ ^ ' 



To estimate the role of the viscosity we need to replace 
(^) with an expression that is relevant at white dwarf 
densities. As a rough approximation, we will use (for a 
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star composed of C ) 



?7 = 9.7 X 10^ 



7/4 



lOSg/cm'' 



lO^K 
T 



1/4 



gem ^s ^ 



(13) 

This expression is deduced from data used by Durisen 
(1973) and estimates the electron contribution to the 
shear viscosity at densities 10'' < p(g/cm ) < 10^. We 
have used the results for non-relativistic electrons since 
they lead to a larger value of 77, in order not to under- 
estimate the role of the viscosity. We have compared 
our approximation to more recent (and detailed) results 
of Itoh et al (1987). The comparison suggests that our 
expression ( p^ ) is a slight overestimate of the strength 
of the electron viscosity (by a factor of less than 5 in 
the temperature/density range of interest). Furthermore, 
we have checked that the electron contribution dominates 
that of the ions in the temperature/density range of in- 
terest. This is true as long as the density is above 10^ 
g/cm^ or so. At lower densities the results of Itoh et al 
(1987) indicate that the ion contribution will dominate 
over the electron contribution. Then we expect (^|) to be 
an underestimate by (again) less than a factor of 5. In con- 
clusion, we think that estimates obtained using ( |l^ ) are 
reasonable, but the potential erring factors of 5 should be 
kept in mind. Such errors would affect the critical rota- 
tion period at the 25% level. Since the main objective of 
the present estimates is to assess whether a more detailed 
study of the role of the r-mode instability in white dwarfs 
is warranted we believe it is reasonable to work at this 
level of accuracy. 

Using (O) we find that 



1.7x10-' 



Mq 

M 



3/4 



R 



O.Oli?, 



17/4 



/ T 



V lO^K 



1/4 



yrs, 

(14) 



and a critical period 



Fc«27 



M 



1/24 



R 



O.Oli?, 



11/8 



/ T 



VlQS K 



1/24 



S . 
. . 

As in the case of neutron stars, the critical period is not 
strongly dependent on the mass. Furthermore, in this case 
it is also weakly dependent on the temperature. The ques- 
tion is: Does the estimated critical period indicate that 
the r-mode instability is relevant also for white dwarfs? 

First of all, we should compare the deduced critical 
period to the mass-shedding limit provided by the Kepler 
frequency. The value of P^, ~ 27 s for a canonical white 
dwarf (at T = 10^ K) should then be compared to the 
Kepler period for a white dwarf: Pk ~ 17 s (that follows 
from fiif = O.Qly/irGp, but if the star is differentially 
rotating at the relevant temperature the limiting period 
can be much smaller than this value). The fact that Pc > 
Pk indicates that the r-mode instability should become 
active before a white dwarf that is spinning up reaches 
the mass shedding limit. 



Moreover, a comparison of our result to the avail- 
able observations of rapidly spinning white dwarfs is 
suggestive. There are no observations of white dwarfs 
rotating faster than the limit set by (|l5|). But the 
fastest spinning white dwarfs, the so called DQ Herculis 
stars (which are magnetized accreting cataclysmic vari- 
ables) (Patterson 1994), have periods quite close to Pc- 
The shortest definite rotation period that has been ob- 
served is the 33 s of AE Aqr. A likely candidate for slightly 
faster rotation is WZ Sge with a period of 28 s. The esti- 
mated surface temperature of these stars is 15000 K, and 
an upper limit (above which the colour of the stars would 
change considerably) is 50000 K (Patterson, private com- 
munication). Just as in the case of neutron stars, the 
interior temperature (which should be used in the r-mode 
estimates) is likely to be considerably higher (maybe two 
orders of magnitude) than the surface value. However, 
since the critical period Pc is very weakly dependent on T 
an increase of the temperature by two orders of magnitude 
only changes the critical period by 20% (at T — 10^ K we 
get Pc ~ 33 s). In Figure |^ we compare our estimated 
critical period to the observations. 

Our estimates suggest that the fastest spinning white 
dwarfs that have been observed are close to the limit 
where the r-modes would be unstable. This is an inter- 
esting result, but it is clear from (|lj) that the r-modes 
would grow very slowly in these stars. To see whether we 
should expect the instability to beable to prevent the star 
from spinning up we compare the growth timescale to the 
estimated spin-up time for a white dwarf accreting at the 
rates observed in the DQ Her. The spin-up timescale can 
be approximated by 



J 
j 



2MR'^n 1 



bVGMR M 



(16) 



which leads to 



1.3 X lO'' 



M 



1/2 



R 



OMR. 



3/2 



30s 



X yrs 



(17) 



Interestingly, the two timescales tg^ and t^u would bal- 
ance at P « 30 s for the observed accretion rate of AE Aqr 
(M = 2 X lO"^°M0/yr). Still, we believe the growth rate 
of the mode is too slow for the instability to be relevant. 
If the star spins up closer to the Kepler period the mode 
grows faster, but at for P = 17 s we still get t^u, « 10^ 
yrs. This means that the r-mode has to be able to grow 
coherently for 100 million years in order for the instabil- 
ity to affect the rotation of a white dwarf, which seems 
unlikely. Hence, we can probably rule out the r-mode in- 
stability from playing any role role in the scenario where 
a weak-field white dwarf is spun up towards the Kepler 
frequency and then undergoes accretion induced collapse. 
As was shown by Narayan and Popham (1989) accretion 
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induced collapse is likely to occur only if the accretion 
rate is higher than 4 x 10~^ Mq /yr. The white dwarf is ex- 
pected to reach critical mass after accreting 0.1 — 0.1 SM©, 
which means that it would collapse after 10*" years or so. 
This would clearly not give the r-modes time to grow by 
the many orders of magnitude required for the instability 
to affect the rotation of the star. 

Finally, a word of caution is in order. Our estimate 
for the viscosity assumed that the stellar fluid was in the 
liquid phase. At low temperatures this may not be the 
case, and the viscosity of the resultant lattice may be 
much higher than (|^). We estimate that the instability 
becomes active unless the true viscosity of white dwarfs 
is more than 400 times larger than (p^). 

5. Final remarks 

We have presented rough estimates that suggest that 
the r-mode instability may be active in accreting neutron 
stars and white dwarfs. Our back-of-the-envelope calcu- 
lations lead to two main conclusions. 

Firstly, the r-mode instability can limit the rotation of 
recycled MSPs provided that they are hotter than 2 x 10^ 
K. We have shown that this suggestion cannot be ruled 
out by present observations for MSP periods and temper- 
atures. We have also compared our predicted critical peri- 
ods to observed data for LMXBs. In particular, we consid- 
ered the systems for which a rotation period in the narrow 
range 3-3.8 ms has been inferred from the beat-frequency 
model for QPOs. This comparison strongly suggests that 
the r-mode instability may be active in these systems, and 
that it can potentially be the agent that limits the rate of 
rotation to the observed periods. Furthermore, we have 
estimated the amplitude of the gravitational waves from 
a neutron star rotating at the r-mode instability limit and 
accreting at the Eddington limit (typical examples could 
be the fastest accreting LMXBs or the so-called Thorne- 
Zytkow objects). The result suggests that this kind of 
source (from within our Galaxy) could be detected by the 
new generation of gravitational-wave detectors. 

There are, of course, many uncertainties associated 
with these suggestions. The most important one concerns 
the basically unknown role of superfluidity. A better un- 
derstanding of neutron star superfluidity and its effect on 
the r-mode instability is urgently needed to firm up the 
present estimates. 

Secondly, we find that the r-mode instability may be 
active in short period white dwarfs. Our order of magni- 
tude estimates (for a C^^ white dwarf of M = Mq and 
R = OMRq) yield a critical period P « 27 - 33 s. This 
number is in interesting agreement with the observed peri- 
ods (> 28 s) of the rapidly spinning DQ Herculis stars. It 
is, however, only an order of magnitude estimate that in- 
dicates that the instability should be taken seriously also 
for white dwarfs and that a more detailed study should 
be attempted. 
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Fig. 1. — The inferred critical period for the r-mode in- 
stability at temperatures relevant to older neutron stars 
(solid lines). The upper line is for a normal fluid star, 
while the lower one is for a superfluid star (only tak- 
ing electron-electron scattering into account, see the main 
text for discussion). The data is for a neutron star with 
M = I.AMq and i? = 10 km. The Kepler limit, which cor- 
responds to P w 0.8 ms for our canonical star, is shown 
as a horizontal dashed line. We compare our theoreti- 
cal result to i) the observed periods and temperatures of 
the most rapidly spinning MSPs (see Reisenegger 1997 
for the data). The surface temperatures are indicated 
as solid vertical lines. The dashed continuation of each 
line indicates the estimated core temperature, ii) the 
observed/inferred periods and temperatures for accreting 
neutron stars in LMXBs, and iii) the recently discovered 
2.49 ms X-ray pulsar SAX J1808.4-3658. 
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Fig. 2. — The inferred critical period for the r-mode in- 
stability at temperatures relevant for white dwarfs (thick 
solid line). The data is for a (Carbon based) white dwarf 
with M = Mq and R = 0.01 Rq. The theoretical result 
is compared to the Kepler limit (thick dashed line) which 
corresponds to P « 17 s, and the observed periods of the 
most rapidly spinning DQ Herculis stars (dashed lines) 
(see Patterson 1994). In order from below, these stars 
are WZ Sge, AE Aqr, V533 Her and DQ Her. The likely 
surface temperature of the fastest rotators is 15000 K , 
and the core temperature (which is the relevant one for 
the r-mode estimates) is likely to be something like two 
orders of magnutide higher. 
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